Absence of static phase separation in the high T c cuprate YBa 2 Cu 3 
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We use 89 Y NMR in YBa2Cu306+;, in order to evaluate with high sensitivity the distribution of hole 
content p in the CuC>2 planes. For y = 1 and y = 0.6, this hole doping distribution is found narrow 
with a full width at half maximum smaller than Ap — 0.025. This rules out any large static phase 
separation between underdoped and optimally doped regions in contrast with the one observed by 
STM in Bi2212 and by NQR in LaSrCuO. This establishes that static electronic phase separation 
is not a generic feature of the cuprates. 



In the cuprates, the hole doping of the Cu02 planes 
induces the unusual features of the phase diagram: high- 
T c superconductivity, strange metal behavior, or pseu- 
dogap. A large body of theoretical work argues that 
this hole doping could be intrinsically strongly inhomo- 
geneous in the planes, forming segregated hole- rich and 
hole-poor regions on a nanoscale [jjj. This phase sepa- 
ration has been proposed to be essential to explain the 
unusual properties of the cuprates, appearing for exam- 
ple as stripes as argued from neutron scattering experi- 
ments |Q . Such proposals have been recently highlighted 
by STM studies which reveal strong inhomogeneities of 
the superconducting properties at the surface of Bi2212 
§ § § § §■ The Berkeley STM group have imaged 
a spatial distribution of the superconducting gap which 
they associated with a distribution of concentration of 
holes ranging from p — 0.1 to p — 0.2 holes/planar unit 
cell®. 

The archetype of the cuprate families is YBa 2 Cu306+ y 
(YBaCuO). If these electronic inhomogeneities do not 
exist in this YBaCuO family, then they result of spe- 
cific disorder in other compounds. However, STM mea- 
surements are not as extensive in YBaCuO due to sur- 
face cleaving problems and oxygen loss in vacuum. The 
existing STM studies display an inhomogeneous or rel- 
atively homogeneous surface depending on the surface 
preparation procedure || Q. These limitations do not 
occur if one uses other local probes such as nuclear mag- 
netic (NMR) or quadrupolar resonance (NQR). The huge 
body of NMR/NQR studies done so far concentrated on 
the p— variation of the average values of specific quan- 
tities (i.e. the static NMR shift K or the relaxation 
time Ti). However, as will be emphasized hereafter, the 
NMR/NQR spectroscopy also allows to determine the 
distribution of these quantities in the bulk samples. NMR 
is then sensitive to the local distribution of electronic 
properties like STM, but not to the specificities of the 
surface. As an appealing example, Singer et al |l0[ ev- 
idenced recently a distribution of T\ over the Cu NQR 
spectrum in bulk LaSrCuO, which can be attributed to 



a distribution of p as large as the one observed on the 
Bi2212 surface. 

We present an NMR static study of the YBaCuO fam- 
ily using 89 Y NMR spectra. This allows us to estimate 
the shape of the hole distribution in the bulk of the com- 
pound. We study two doping compositions (underdoped 
06.6 and slightly overdoped O7). At O7 the oxygen chains 
reservoirs are full, thus well ordered. At 06.6, the various 
properties are weakly oxygen content dependent. So both 
compositions should produce a minimal disorder of dop- 
ing in the planes. We show indeed that our results reveal 
a very narrow hole distribution. These results will be 
compared with other experiments in the various cuprate 
families. 

We synthesized a large sample batch of single crys- 
tal grains of YBaCuO. It was oxidized at T = 350°C, 
and then slowly cooled under oxygen atmosphere from 
350°C to 270°C at a rate of OAK /hi. After this proce- 
dure, the sample has a nominal maximum oxygen con- 
tent y = 1 with a reduced T c — 89.1 K with respect 
to optimal doping (T c — 92. 5K). A part of this batch 
was reduced to y = 0.60 ± 0.02 (T c = 55.8K) under pri- 
mary vacuum (P = 0.13 mbar) at T = 385°C . The 
oxygen content was measured by thermogravimetry, i.e. 
measurement of weight loss during the deoxidisation pro- 
cess until equilibrium is reached in the above conditions. 
For both samples, the crystallites were then aligned in 
stycast epoxy in the field H ext (~ 7.5 Tesla) of the NMR 
spectrometer pT| . This field corresponds to a reference 
NMR resonance 89 v = 15634.67 kHz for a liquid YC1 3 so- 
lution. Spectra were obtained for H ex t perpendicular to 
the c cristallographic axis using a standard tt/2 — tt pulse 
sequence and Fourier transform of half of the spin echo 
Jl2[ | . The repetition time of the pulse sequence was long 
enough at each temperature to recover the saturation sig- 
nal for any value of p. For a given T, the maximum of 
Ti(p) varies from 17 to 110 sec when decreasing T from 
300K down to 120-K" Q. 

89 Y NMR probes sensitively any doping distribution in 
the Cu02 planes because the Y nucleus is coupled to Cu 



1 



of its adjacent CuC>2 planes through the oxygen orbitals 
|l3f . This results in a shift of the NMR line related to the 
Cu spin susceptibility x given by 



J A(T) = 



89 



Ms 



(1) 



Here, 89 A(T) is strongly dependent on hole content 
because of the strong hole doping dependence of x(T)- 
In contrast, the T— independent chemical shift 89 S and 
hyperfine coupling Ahf = —1.95 kOe/fj,B between the 
89 y and one Cu do not change much with hole doping 
(hb is the Bohr magneton) |13j. This relation is illus- 
trated in Fig.l which displays the shift of the peak value 
of the 89 Y NMR line. It has been taken for a series 
of samples prepared from a single batch for different y. 
The shift 89 A (T) displays the well known behaviour for 
x(T). It is nearly constant above optimal doping. It dis- 
plays a marked decrease at low T which is the signature 
of the pseudogap for the underdoped cases. The oxy- 
gen content y can be converted into the hole content per 
plane using the parabolic dependence of T c jl4| given by 
p = 0.16 ± -T c (y)/T™ ax )/82.6. The corresponding 
p and experimental T c are plotted in the inset of fig.l. 
Note that near optimal doping, the large p dependence 
of 89 A makes it a sensitive probe to any doping variation: 
for example, at T = 120 K, a change in p as small as 0.01 
corresponds to a change in 89 A" of about 20 ppm, easily 
detectable. In the underdoped regime, the sensitivity to 
the doping is smaller but still sizeable. 

The NMR spectrum of a given sample consists of a his- 
togram of the NMR shifts throughout the sample. Note 
that, as STM, the NMR local probe is sensitive to vari- 
ations nearly on the atomic scale, as each Y nucleus is 
coupled only to its 8 nearest neighbour Cu sites. Then, 
from the correspondence shift-doping of fig.l, any local 
distribution of doping will lead to a similar distribution 
in the shift, hence to a broadening of the spectrum. 

Let us now present the experimental 89 Y NMR spectra 
for y — 0.6 and 1 from which we will extract quantita- 
tively the actual doping distribution. The spectra are 
plotted for T = 120 K and 300 K for the optimally and 
underdoped samples in fig. 2. As can be observed imme- 
diately there is no significant overlap between the two 
spectra at T = 120 K . This already proves that the 
O7 (p = 0.18) sample contains no appreciable amount 
of underdoped nanoscale regions (equivalent to p = 0.1) 
and vice versa. In fig. 2 we also display the hole doping 
scales deduced from the relation between p and 89 A ob- 
tained from fig.l. From these hole doping scales it can 
be seen roughly that the full width at half maximum 
FWHM of the doping distribution is necessarily smaller 
than Ap ~ 0.04 at optimal doping and 0.02 for the under- 
doped case. This clearly proves without any further anal- 
ysis that the observed doping distribution in YBaCuO 
is much smaller than in Bi2212 and LaSrCu30 where 
Ap ~ 0.1. 

For a more quantitative analysis, we need to consider 



besides the distribution of x(T) any distribution of the 
chemical shift 6 and hyperfine coupling Ahf which en- 
ter Eq.|l]. Such distributions will come from any local 
structural disorder. These distributions are fully respon- 
sible for the FWHM AK = 32 ± 5 ppm of the undoped 
compound for T > Tn 15 1. Indeed, for oxygen content 



yC 0.15, no hole distribution broadening is expected 
. This results from the fact that xoxygens introduced 
in the Cu reservoir layer mainly convert the two adjacent 
3d 10 Cu(l) into 3d 9 and do not yield any hole doping of 
the CuC>2 planes, so that p is strictly 0. It is confirmed 
by the fact that the spectrum has the narrowest width 
amoung all dopings and is not found significantly depen- 
dent on oxygen content up to y = 0.15. 

For higher y, if AS and AA hf are the FWHM of the 
respective distributions assumed to be uncorrelated and 
gaussian for simplicity, this will lead to an additional 
broadening with FWHM AK(y, T): 



AK 2 (y,T) = AS 2 (y) 



AA^(y,T) 



(2) 



where AK hf = AA hf x(T). We note that the y = 
compound is the more ordered as it is tetragonal with no 
twin boundaries and empty chains, hence AS and AAhf 
should increase at higher y. As these broadenings add 
to the doping induced broadening, taking their y = 
estimates at higher dopings will lead to an overestimate 
of the doping distribution. As AKhf is T-dependent, 
we can evaluate AAhf and AS separately. Indeed, at 
y = 0.6, x( T ) doubles between 120 K and 300 K so that 
AKh f doubles as well, whereas the FWHM increases only 
by 11% as seen in fig. 2. As this T-dependence might as 
well be due to the doping distribution, this leads to an 
upper bound AAhf /Ahf ^ 0.13. Such an upper bound 
would be explained in a naive hybridization computa- 
tion by a random displacement of Cu by 0.04 A from 
its ideal position, a typical value from Rietveld measure- 
ments. We will consider the two extreme cases: 

(i) AAhf /Ahf = corresponding to AS = 32 ppm 

(ii) AAhf /Ahf = 0.13 corresponding to AS = 29 ppm. 
In order to extract the actual doping distribution, we 

start from a distribution of oxygen doping P(y). We 
convert it for different temperatures into a shift distri- 
bution through the phenomenological 89 A vs y variation 
deduced from Fig.l. We convolute this A distribution 
with the 5 and Ahf gaussian distributions with FWHM 
either (i) or (ii). This process is iterated until an opti- 
mal P(y) is found to fit the experimental spectra for all 
temperatures with no additional free parameter. 

In the underdoped y — 0.6 compound, this fitting pro- 
cedure is limited by the fact that the NMR shift is nearly 
insensitive to the doping distribution for y < 0.6 as seen 
in fig.l. However, the thermogravimetry procedure dur- 
ing deoxydation constrains the measure of the average 
oxygen level y = 0.60 ± 0.02. We thus assume P{y) to 
be symetric around y. The best fit is then obtained for a 
gaussian distribution P(y) = exp (—(y — 0.6) 2 /er 2 ) with 
a = 0.05 for (i) and a = 0.1 for (ii). This oxygen dis- 
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tribution P(y) and the corresponding hole distribution 
V(p) are plotted in fig. 3 together with the usual T c di- 
agram. The corresponding simulations fit perfectly the 
experimental ones at all temperatures (examples of fits 
are given as dotted lines in fig. 2). 

For YBaCuC>7, the chains are completely full. Further 
oxidation is prohibited. Therefore the distribution P(y) 
is expected to be non symmetric fl7|| . To take this lim- 
itation into account, we model P(y) as a convolution of 
a gaussian by exp(— \y — 1| /A±) where A± are allowed 
to differ for y > 1 and y < 1. The best fit is found 
for a = 0.01, A_ = 0.06, and A+ = 0.025 (i) or 0.004 
(ii). P(y) and the corresponding fit of the spectra are 
respectively plotted in fig. 3 and fig. 2. The high asym- 
metric shape towards y < 1 found for P(y) confirms that 
no large overdoping of the planes is produced locally, as 
expected. In order to illustrate the high accuracy of our 
method, we choose P{y) slightly broader than our best 
fit, plotted as the dotted line in the upper panel of fig. 3. 
The corresponding simulation plotted as a dashed line in 
the upper panel of fig. 2 clearly fails to fit the experimen- 
tal spectrum. 

In summary, our results demonstrate that the maxi- 
mum possible distribution of doping is quite sharp with 
typical width Ap ^ 0.025 for optimal doping and Ap ^ 
0.01 for the underdoped sample. The distribution is much 
smaller than in the LaSrCuO and Bi2212 cuprates. In 
LaSrCuO, at optimal doping, Singer et al. find Ap ~ 0.09 
as figured by the arrow in fig. 3 En]. Another NMR study 
||18t using 17 and 63 Cu NMR allowed to evidence a short 
length scale spatial modulation in underdoped LaSrCuO 
compatible with Ref. §j(J. In Bi2212, the STM experi- 
ments reveal regions on the surface with the usual super- 
conducting gap while others display a pseudogap || [Q ||] 
H |U . This has been interpreted to be due either to some 
disorder effect Q or to superconducting and insulating- 
like underdoped regions [f| @. The corresponding 
doping distribution of width Ap — 0.085 extracted in 
Ref. |5j is plotted in fig. 3. Similar widths are found for 
an optimally and an underdoped compound in Ref. [Q. 
Therefore, YBaCuO appears much more homogeneous 
than both LaSrCuO and Bi2212. In other cuprates, no 
key experiment was performed to probe locally the dop- 
ing distribution, to our knowledge. However, in Hgl201 
and T12201, typical 17 NMR widths are similar to those 
in YBaCuO and much smaller than in the Bi and La 



compounds |19|. Following the analysis done above and 
the proportionality between 17 and 89 Y shifts, we then 
expect the Hg and Tl family not to exhibit any strong 
doping distribution either. 

In the LaSrCuO family, the large distribution seen by 
NQR occurs in the bulk of the material and might be 
associated with tilt of the oxygen octaedra, buckling of 
the planes, or stripes. In the Bi material, the STM re- 
sults might be only a surface specificity. In YBaCuO, the 
doping distribution P(y) measured here is not only much 
narrower but might even have only a macroscopic ori- 
gin. As we used powders of cristallites with sizes smaller 



than 30 /xm, any oxygen gradient within each cristallite, 
or a correlation between y and the cristallite size could 
lead to the observed P(y). At O7, such oxygen gradients 
are explained by the fact that oxygen diffusion between 
chains becomes limited at low temperature. This natu- 
rally leads to the asymmetry of P(y) towards y < 1 seen 
in fig. 3. This effect has been observed systematically in 
the various studies published by our group using 89 Y, or 
17 NMR spectra in many different YBCO7 powders. At 
06.6, the observed broad P(y) is very narrow when plot- 
ted versus p (fig. 3). Indeed a change in oxygen content 
does not strongly modify the actual hole doping of the 
planes near y — 0.6 hence leading to the plateau of T c . 
This is to be associated with the fact that extra oxygens 
at such composition occupy empty chains, and do not 
modify the hole content, in analogy with the situation 
at y = 0. Therefore, the actual distributions would be 
much narrower on a submicron size region similar to the 
one sampled by STM. For both dopings, the key factor 
to the good homogeneity encountered is then probably 
the presence of the chains. The existing chain disorder 
does not lead to a sizeable distribution of hole content in 
the planes. Furthermore, those chains are probably suffi- 
ciently far from the planes not to stabilize a charge seg- 
regation. This specificity makes YBa 2 Cu306.6 or 7 some 
of the best prototypes of clean homogeneous cuprates. 

In conclusion, the nanoscale static segregations ob- 
served so far cannot be considered as an intrinsic phe- 
nomenon common to all cuprates. Contrarily to the ob- 
servations in Bi2212 and LaSrCuO, the maximum hole 
doping distribution Ap found in YBCO is small enough 
to conclude that our YBCO samples do not consist in 
interleaved regions with qualitatively different physical 
properties (metal versus insulator, with or without pseu- 
dogap, etc). The only remaining possibility for charge 
segregation or stripe-like scenarii to apply inYBaCuO 
is then a dynamical process where the phase separation 
would have a lifetime smaller than our timescale of obser- 
vation. For the present experiment, it corresponds to the 
inverse spectral width, typically 1 msec. Such dynamics 
might have been detected in inelastic neutron scattering 
experiments, but the existing results are still under de- 
bate |o| gij. 

Fig. Legend : 

fig.l : 89 Y NMR shift S9 K for H ext ±c for different 
oxygen contents y versus temperature. In the inset, T c 
for each sample is plotted versus hole doping with the 
same symbol as the one used for K. 

fig. 2 : 89 Y NMR Fourier Transform spectra are plotted 
at different temperatures and dopings in YBa2Cu306+y 
for H ex t-Lc (black solid lines) with arbitrary normaliza- 
tion. In each panel, a doping scale gives the relation 
between the shift 89 K and hole doping p. Dotted lines 
are simulated spectra (see text). At Oj:T = 120 K, an 
additional dashed line is plotted, which corresponds to 
the simulation made using the dotted distribution of the 
upper panel of fig. 3. 
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fig. 3 : The full black dots represent T c versus oxy- 
gen content (upper panel) and hole doping (lower panel) 
for YBa 2 Cu306+ a . Distributions of oxygen content P(y) 
and the corresponding distribution of hole content V{p) 
were used to fit the spectra for y = 1 and y = 0.6 sam- 
ples (see text). The dark and light gray enveloppes are 
obtained with a minimal distribution of chemical shift 
and hyperfme coupling within assumptions [i) and (ii) re- 
spectively. The dashed distribution and the arrow shown 
in the lower panel represent the distributions measured 
in Bi2212 § and LaSrCuO Jl(J. 

*present adress : Dept of Physics, IIT Bombay 400076 
India. 
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